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Abstract 

We address the question of whether atomic bound states begin to sta- 
biUze in the short ultra- intense field limit. We provide a general theory 
of ionization probability and investigate its gauge invariance. For a wide 
range of potentials we find an upper and lower bound by non-perturbative 
methods, which clearly exclude the possibility that the ultra intense field 
might have a stabilizing effect on the atom. For short pulses we find almost 
complete ionization as the field strength increases. 
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1 Introduction 



Fermi's Golden Rule, as one of the central elements in quantum mechanics, has 
served for many years for the understanding of photoionization rates of atoms 
in weak radiation fields. Its origin is, however, perturbative and therefore when 
applying very intense fields (with intensities which are greater or of the order 
of one atomic unit 3.5 x lO^^W/cm^) one leaves its range of validity. With the 
advance of laser technology this high intensity region has become accessible to 
real experiments in the form of laser pulses of Ips or less, at frequencies ranging 
from the infrared to the ultraviolet [||]. The prediction of atomic ionization rates 
are of practical importance for instance in the study of gas breakdown 0. 

In order to treat the new regime, several alternative approximation methods 
have been proposed. On one side P, ^, Q they are based on a perturbation 
around the Gordon- Volkov solution |^ of the Schrodinger equation. The question 
of convergence of these series and their precise range of validity has not yet been 
put on firm grounds. Despite these problems, these methods have been applied to 
find numerical solutions for the ionization probabilities. On the other side there 
exist a vast number of numerical studies, which make use of numerical solutions 
of the Schrodinger equation, high-frequency approximations or the Floquet 
approximation H]. Most computations have been carried out in one dimension 



T0| , in the hope that the essentials of the full three-dimensional physics are 



already present in this simplified situation. There exist arguments which put 
them in question |jTT], |T^, since in comparison with the full three dimensional 
situation, they do not account for the full angular dependence and may provide 
misleading results. Recently there have been full three dimensional computations 
||T3| , p!5| , p!2| , |T6|, |1^] . But the complete problem has not been solved yet and as 



it is pointed out in [|18[, "even the simplest one-electron atom in an intense laser 
field presents too great a challenge for truly ab inito numerical work, and a variety 
of compromises have been developed" . These compromises are partly located in 
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the numerical methods themselves, but put partly constraints onto the physics, 
such as for instance the introduction of mask functions or the approximation of 
the continuum by cutting of the high energetic states. 

Several authors claim to have found the very surprising and counter-intuitive 
result, that the bound state of the atom stabilizes as the field strength increases 



Tj, 0, |T7|. Similar results have also been obtained by many authors for 



the one dimensional situation 0. In fact these findings are so surprising that "a 
dramatic shift in viewpoint is required to explain the physics of atoms in very 



strong laser fields " [|1^]. We shall comment on these results below and for the 
moment refer the reader to the review article on these findings by Eberly and 
Kulander [|Ty| and one by Geltman |]2D[, who takes the opposite point of view 
that atomic bound states do not stabilize as the field strength is increased and 
who asserts that the "conventional interpretation of the theory of the interaction 
of radiation and atoms is quite sound even in this regime". The latter point of 



view is also supported in ||21 



Evidence for atomic-stabilization in superintense laser fields has also been 
obtained from the study of several classical dynamical systems p2| |. 

Up to now there exist no data for intensities of one atomic unit, i.e. in the 
high intensity region for which the theoretical predictions are made, such that 
the controversy could be settled from the experimental side. So far there exist 
some experiments for lower intensities lO^^W/cm'^, which provide evidence for 



some sort of stabilization |23|. 

The controversy is mainly based on numerical results and a detailed theoretical 
analysis of the problem which involves analytic expressions only does not exist so 
far. The main intention of our paper is to provide an alternative approach to the 
matter. We consider the Schrodinger equation for an atom in a linear polarized 
electric field, 

^^ = (-A/2 + V + z■E{t))^p, 
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where E{t) stands for the intensity of the field and is supposed to have finite 
duration (for instance, 1 ps = 4.17x10'^ a.u.). We do not specify E{t) in more 
detaih it can be, for instance, a pulse which contains a number of optical periods 
(the frequency of which is determined by the frequency of the laser) possibly with 
some turn-on and turn-off parts. Such kind of pulses were used in the search for 
possible suppression of ionization. We note that 1 ps pulses have a duration 
comparable with a classical Kepler period for the highly excited Rydberg states 
{Tk = 27m^). Another example are half-cycle pulses with duration of about 500 
fs, generated in the experiments of Jones et al. ^4 . However, the maximal 



intensity reached in these experiments was about 10"^ a.u., so that these pulses 
are ultimately far from the ultra- intense limit. 

We suppose furthermore that the wave function ip{x, t) is given by the bound 
state wave function of —A/2 + V before the pulse is turned on. One can easily 
estimate that relativistic effects might be appreciable as soon as E{t) is so strong 
that the classical theory predicts electron velocities approaching the speed of 
light, or more precisely when in atomic units the electric field strength time the 
frequency is of the order of the fine structure constant. According to the estimates 
in |jl8| this occurs for typical frequencies for laser intensities E ^ lO^^W/cm^. 



Our results below show that atoms do not become resistant to ionization when 
exposed to short ultra-intense laser pulses. Our statements are of qualitative na- 
ture, in the sense that they provide upper and lower bounds and do not predict 
precise values of the ionization probabilities. The methods we use cover all pos- 
sible pulses, i.e. also those which are very popular in the literature with smoooth 
turn on and off. Our arguments cover all frequency regimes, including the high 
frequencies for which stabilization is supposed to occur. For pulses which are not 
switched on smoothly, our results typically hold for very short times of the order 
of one atomic unit. With a smooth switch on of the pulse one may extend the 
region of validity. We provide expressions for two upper bounds, (3.2) and (3.14) 
valid in the region when {J^ E{t)dtf' /2 (the classical energy transfer of the pulse) 



is smaller than the binding energy and the other valid without restriction. The 
lower bound holds when {£ E{t)dtf /2 > -E. 

The paper is organized as follows. In Section 2 we formulate a general theory 
of ionization probability and prove its gauge invariance. We also make con- 
tact to the various approximation methods based on perturbative expansions. 
In Section 3 we briefly discuss these methods in the context of quantum me- 
chanical one-particle Stark Hamiltonians and provide the proofs for the upper 
and lower bounds for the ionization probability for a wide range of one-particle 
potentials, which in particular include all potentials appearing in atomic and 
molecular physics. In Section 4 we state our conclusions. In Appendix A we pro- 
vide an upper bound for the Coulomb potential and in Appendix B we optimize 
this bound for the ground state of the hydrogen atom. 

2 The general theory of ionization probabihty 
and its gauge invariance 

In this section we will give a general discussion of the ionization probability 
and its gauge invariance. Gauge invariance is of course necessary for observable 
quantities and it is conventional wisdom for the case of the ionization probability. 
However, we were not able to locate an explicit reference with a proof and will 
therefore include a discussion on this issue. We will relate our arguments to 
familiar concepts in scattering theory and explicitly discuss its relevance in the 
context of the Stark Hamiltonian. In the last part of this section we will show 
the gauge covariance of time- dependent perturbation theory. In order to convey 
the general ideas we will avoid bulky mathematical notations in this section. 

Let H{t){—oo < t < oo) be a general time dependent selfadjoint Hamilto- 
nian in some Hilbert space H and let U {t, t') denote the resulting time evolution 
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operator from t' to i.e. U{t,t') satisfies 

idtU{t,t') = H{t)U{t,t') 
U{t,t')U{t',t") = U{t,t") (2.1) 
U{t,t) = 1 

for all t, t' , t". In the context of the Stark Hamiltonians, U{t, t') exists for all t, t' 
and is unitary (see below). 

Assume now that H{t) approaches an operator for t ^ oo and i7_ for 
t —>■ — oo, i.e. 

H+ = lim H(t) 

H_ = lim H{t) (2.2) 

t— >— oo 

holds in a suitable sense. It is important to note that we do not assume if+ to 
equal H_. In fact, for the Stark Hamiltonian in certain gauges, these operators 
will in general differ (see below). In analogy to the scattering matrix (see below) 
we define the abstract S-Matrix to be the following weak limit, i.e. the limit for 
matrix elements (if it exists) 

S = ^lim^ exp itH+ ■ U{t, t') ■ exp -it'H_. (2.3) 

t'->-oa 

In particular S exists trivially and is unitary if H(t) becomes stationary for all 
large |t|, i.e. if H[t) = for all t > t+ and H{t) = H for all t <t^ for suitable 
finite In this case we will call H{t) a finitely pulsed Hamiltonian. S then 

takes the form 

S = exp itH+ ■ U{t, t') ■ exp -it'H_ (2.4) 

for all t > t+ and all t' <t^. In particular S is then unitary. 

Let P+ be the orthogonal projection onto the subspace spanned by the bound 
states of H+. P_ is defined analogously in terms of Then for any normalized 
state ip in the range of P_ its ionization probability is defined to be 

m^Wil- P^)Si;r. (2.5) 
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Here denotes the Hilbert space norm, i.e. WipW^ = For the case when 



= this agrees with the definition used in ^ . For a finitely pulsed 
Hamiltonian we have 

m = \\{l-P^)U{U,t.W (2.6) 

whenever is a. bound state of H_. 

Abstract gauge transformations are now introduced as follows. Let A{t) {—oo 
<t< oo) be a one parameter family of unitary operators (suitably differentiable 
in t). If il){t) is a solution of the Schrodinger equation 

then ip'{t) = A{t)tp(t) is a solution of the equation 

idt^P'it) = i{dtA{t))^it)+A{t)idt^{t) 

= tidtAit))Ait)-'^p'{t)+A{t)H{t)A{t)-'^p'{t) (2.7) 

with 

H'{t) = i{dtA{t))A{t)-^ + A{t)H{t)A{t)-^ (2.8) 

being formally selfadjoint. If U'{t, t') is the time evolution operator for H'[t) then 
obviously 

U\t, t') = A{t)U{t, t')A{t')-\ (2.9) 

We note that the set of all gauge transformations form a non- commutative group 
under the obvious multiplication rule {AiA2){t) = Ai{t)A2{t), with unit l{t) = 1 
and inverse A~^{t) = A{t)^^. The familiar interaction picture used in scattering 
theory is now a special case. Indeed, to be more specific assume H{t) to be of 
the form H{t) = Hq + Hj{t), where Hq is the "free" Hamiltonian and Hi{t) the 
(possibly time dependent) interaction Hamiltonian. Set 

A{t) =expitHo. (2.10) 



Then 

H'{t) = H',{t) (2.11) 
is the Hamiltonian in the interaction picture with 

Hj{t) = exp itHo ■ Hi(t) ■ exp -UHq. (2.12) 

In this case the hmit (if it exists) 

S(H,Ho)^ hm U'(t,t')= hm ex.p itHo ■ U(t,t') ■ ex.p -it' Ho (2.13) 

t— » + oo t— > + oo 
t'->-oo t'^-oo 

is called the scattering matrix (S-matrix) for the pair {H,Ho). Let us elaborate 
briefly in what sense there is an analogy between the quantity (2.3) and the 
case of ordinary potential scattering, i.e. where Ho is —A/2 and Hi{t) = V 
is a potential. In the latter case H — Ho + V is compared with Ho in spatial 
regions far out, i.e. where V is small and where the incoming and outgoing wave 
packets are located for large \t\. In the case we are interested in presently, H{t) 
is compared with H^ for large positive times and with H_ for large negative 
times. To summarize: In one case one compares Hamiltonians for large spatial 
coordinates and in the other case for large time coordinates. 

To apply this general concept of gauge transformations and gauge covariance 
to our discussion of ionization, assume now in addition that A{t) approaches 
suitably unitary operators and A_ when t — > +oo and t — > — oo, respectively. 

Then H'{t) (see (2.8)) approaches 

H'^^A+H+A-^ (2.14) 

and 

H'_ - A_H_AZ^ (2.15) 

as f — > +00 and i — > — oo respectively. 

By formal manipulations we therefore have 

S' = A+SAZ\ (2.16) 
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In particular if H{t) is finitely pulsed and if in addition A{t) is stationary for all 
large then H'(t) is also finitely pulsed, S' exists, is unitary and (2.16) holds. 

In general, by (2.14) and (2.15), if P4 the orthogonal projections onto the 
space spanned by the bound states of H'j_, we have 

P± = A±P±Al\ (2.17) 

In particular A_tlj is in the range of PL if t/j is in the range of P_. Inserting (2.16) 
and (2.17) gives the desired gauge invariance in the form 

I'(A_^|;) = \\(l-P:,)S'A_^j;r 

- \\A+{1- P+)A-^A+SAz^A_ijf (2.18) 

= \\{i-p^)s^r^m, 

since by assumption A^ is unitary. In the example we will be interested in 
A_ — 1. In this case H'_ — H_ and (2.18) takes the simpler form /'(V') = -^(V')- 
Prom the proof we see that gauge invariance is an important regulating principle 
in the following sense. One has to choose the projection P+ in (2.5) and not P_ 
in order to obtain gauge invariance. 

We now apply these concepts to the theory of the time dependent Stark 
Hamiltonian. In order to make notations more transparent we choose a linearly 
polarized electric field, which does however not limitate our discussion since more 
general fields may be simply obtained by replacing z ^ x and E{t) — > E{t) 
{x E Ei^jE E Si^), such that particular other choices, like for instance, circular 
polarized light may be easily be derived from there. We do, however, assume a 
dipole approximation, such that the electric field becomes a function only of time 
and thus is independent of space. Then using atomic units h — e — rrie — c- a — 1 
we consider on the Hilbert space L'^{IB?, d^x) the 3 time-dependent Hamiltonians 

H^{t) = -^ + V + z-E{t) 

H2{t) = \{-iV-h{t)e,f + V (2.19) 
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Hsit) = -^ + V{x-c{t)e,) . 

Here V is an arbitrary potential, is the unit vector in the 2;-direction and 
V{x—y) is the shifted potential, i.e. the multiplication operator on wave functions 
given as (y{x — y)ip){x) = V{x — yjip^x). Also E(t) is the electric field, assumed 
to vanish unless < t < r (i.e. t_ = and t+ = r > in the notation above). 
Apart from this condition the pulse E{t) may be arbitrary. We only make the 
mathematical restriction, that E{t) is piecewise continuous, which means that 
the pulse may have jumps and all commonly used enveloping shapes, for instance 
cosine squared, smooth adiabatic turn on and off, etc., are included. Then the 
following quantities a{t),b{t) and c{t) are well defined 

b{t) = f E{s)ds (2.20) 
Jo 

c{t) = [ b{s)ds = tb{t) - / sE{s)ds (2.21) 
JO Jo 

a{t) = - fb{sfds . (2.22) 
2 JO 

Note that b{T)ez describes the classical momentum transfer of the pulse, such 
that |&(t)^ is the classical energy transfer. Also c{t)cz is the classical displace- 
ment caused by the pulse. Then H2(t) is obtained from Hi{t) by the gauge 
transformation 

y42<_i(t) = exp ib{t)z. 



Hiit) is obtained from H^{t) by the Kramers-Henneberger transformation [^, 

Ai^-iit) = T(t) = exp —ia(t) ■ exp —ib(t)z ■ expic{t)pz. (2.23) 

Therefore we call Hs{t) the Hamiltonian in the Kramers-Henneberger gauge. 
We note that a corresponding transformation in Quantum Electrodynamics was 
introduced by Pauli and Fierz already in 1938 but with a different motivation 



|31|. The Hamiltonian Hi{t) is usually referred to as the Hamiltonian in the 
length- or electric field gauge, whereas H2{t) is denoted as the Hamiltonian in 
the velocity-, radiation- or Coulomb gauge. 



As a consequence H2{t) is obtained from H^lt) by the gauge transformation 

^2^3^ = A2^iit)Ai^3{t) = exp -ia{t) ■ exptc{t)p,. (2.24) 

The general gauge transformation Aj^i{t) for Hi{t) Hj(t) is then obtained 
from the rules 

Aj^i{t) = Ai^j(t) 

A.^kit) = Aj^i{t)Ai^kit). (2.25) 

Hi{t) and H2{t) are finitely pulsed. Note however that in general H^lt) is not 
finitely pulsed and hence has not always a proper limit as t +oo. This is 
due to the fact that b{t) is constant for t > r such that c{t) grows linearly in 
t > T whenever 6(r) ^ 0, as is apparent from ( [^.21 ). Nevertheless the Kramers- 



Henneberger gauge is quite useful as we shall see below. These observations are 
related to the fact that A2^i{t) becomes stationary for large t but in general not 
A2^3{t) and Ai^^lt). Note that by assumption on E{t), Aj^i{—oo) = 1 and thus 

= = H2,- = H^,- =Ho + V. (2.26) 

However, in general -f^2,+ 7^ -^2,-- 

The case = is of special interest, since it corresponds to the situation in 
which the Schrodinger equation admits an exact solution, which is usually referred 
to as the Gordon- Volkov solution 0. Call the resulting operators HQi{t){i = 
1, 2, 3), such that in particular ifo,3(^) = —A/2. The kernels of the resulting time 
evolution operators Uo^i(t,t') can be calculated explicitly. Indeed, we start from 
the familiar relation for the free particle evolution operator (see e.g. [^) 

Uo,3{x,t; X ,t') = {x\Uo^3{t,t')\x) = {x\ expi{t — t') — \x) (2.27) 

expi^ (2.28) 



{27ii{t - t')f/^ 2{t-t') 
Obviously 

(f I exp ih{t)z\x) = exp ib{t)z ■ S'^^x — x') 
(x| exp ia{t)\x') = exp ia{t) ■ S^{x — x) 
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and 

{x\expic{t)pz\x') = {x — c{t)ez,x') = S'^i^x — x' — c{t)ez) 

such that 



{x\T{t)\x') = (f'|T(t)-i|f) 

= exp-ia{t) ■ exp-ib{t)z ■ 5^{x- x' - c{t)ez). (2.29) 

By (2.9) we immediately obtain the following (again well known) relations (see 
e.g. il,||,|,|ll). 

f/o,i(f,t;f' ,t') = (27r2(t - t'))y^ expi{a{t') - a{t)) 

./,//N ,/x /N (x — c(t)ez — x' + c(t')ezY 



Uo,2{x,t;x,t') = (-27ri(t - t'))^/^ expi{a{t') - a{t)) 



X exp t 



2{t- t') 
{x-c{t)ez -x' + c{t')ezf 



2{t- t') 

The kernel of f/o,i(t, t') is often called the Gordon- Volkov propagator (see e.g. 
I, ^,[T1). 

We finally give a discussion of time dependent perturbation theory and its 
gauge covariance. Returning to the general set-up, let H(t) be a " perturbation" 
of K{t). If W{t,t') denotes the time evolution for K{t), we have the generalized 
Du Hamel's formula (see e.g. [Q) in the form 



U{t,t') = W{t,t') - [' ^[U{t,s)W{s,t')]ds 



ds 



W{t,t')-if U{t,s)[H{s) - K{s)]W{s,t')ds. (2.30) 
Jt' 



We recall that in terms of the sometimes more familiar Green's function 

GH{t,t') = -iU{t,t')6{t-t') , 

which satisfies 

{idt-H{t))GH{t,t') = 5{t~t')^ , 
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relation (2.30) takes the form 

+00 

Gnit, t') = GK{t, t') + j Gnit, s) [H{s) - K{s)\ Gk{s, t')ds. (2.31) 

—00 

Similarly one derives the relation 

U(t, t') = W(t, t') - I t Wit, s) [H{s) - K{s)] U(s, t')ds (2.30') 
and hence 

+00 

Guit, t') = GKit, t') + j GK{t, s) [H{s) - K{s)] Gnis, t')ds. (2.31') 

—00 

We recall that Du Hamel's formula in the form (2.31) and (2.31') is the time 
dependent version of the Lippmann-Schwinger equation (see e.g.[^ ^) in the 
case where both H{s) and K{s) are actually time-independent. Indeed, the 
Lippmann-Schwinger equation may be obtained from these relations by taking 
Laplace transforms. 

The equations (2.30), (2.30'), (2.31) and (2.31') may be iterated by intro- 
ducing the left hand side into the right hand side, resulting in a "power series 
expansion" of U{t,t') in powers of H{s) — K{s) and involving W of the form 

00 

U{t,t') = J2Un{t,t') (2.32) 

n=0 

with Uo(t, t') = W{t, t') and with an analogous expansion for Gh- Let us consider 
what happens under gauge transformations. Du Hamel's formula is compatible 
with gauge transformations in the sense that the relation 

f/'(t, t') = W'{t, t') - i f U'it, s) \H'{s) - K'{s)] W'is, t')ds (2.33) 

Jf 

either follows directly for the pair H'{t), K'{t) or by applying the gauge transfor- 
mation to (2.30) and using (2.8) and (2.9). This implies in particular that 

U:,{t,t') = A{t)Un{t,t')A{t')-' (2.34) 



12 



for all n. This applies in particular to the choices H{t) = Hi{t) and K(t) = Ho^i{t) 
given by (2.19) and the gauge transformations which relate them. However, 
some of the approximation methods used for high intensities, on which we shall 
comment more below, use the fact that one can decompose the Stark Hamiltonian 
in two different ways, that is either treating the potential or the term related to 
the electric field as a perturbation. Hence one obtains two versions for (2.30), 
which one may combine iteratively. The series generated in this manner in general 
does not respect gauge invariance order by order. A discussion of this problem 
and a remedy for restoring gauge invariance by including some terms of next 
order, thus leading to cancelations, may be found in |3^, |3^ . 

In the case of the interaction picture (see above) and with the choice K{t) = 
Ho such that K'{t) = 0, the iteration of (2.33) in powers of H'j{t) = H'{t)-K'{t) is 
just the famous Dyson series of the S-Matrix S{H^ Hq) in the limit t —* +oo, t' 
— oo. In the time independent context of the Lippmann-Schwinger equation (see 
above) this corresponds to iterating the Lippmann-Schwinger equation to obtain 
the Born series for the S-Matrix (see e.g. [0]). In the context we are presently 
interested in, such series expansions for the time evolution operator of finitely 
pulsed Hamiltonians lead to a series expansion for the ionization probability when 
inserted in (2.6). In the next section we will discuss the various approaches used 
so far for the Hamiltonian Hi{t) and its gauge transforms H2{t) and H^lt) given 
in (2.19). 



3 Ionization of atoms in strong, short electric 
fields 

Using the notions of the previous section we start with a review and comparison 
of methods and results obtained by previous authors. Then we relate this in a 
first step to a new, rigorous upper bound on the ionization probability, valid for 
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all small r and small classical momentum transfer 6(r) and small displacement 
c(r) (see the upper bound below). This result is also compared with another 



rigorous upper bound previously obtained by two of the authors (V.K., R.S.) [|| 
as well as with a result obtained using a time-energy uncertainty relation given 
by Pfeifer . Secondly we prove a lower bound below, valid for all small r and 
all large &(t), which in particular proves the absence of stabilization. 

Taking H{t) = Ho + V + z-E{t) = Hi{t) and K{t) = Ho + V in (2.30), the re- 
sulting perturbation series in the time dependent interaction H(t) — K{t) = zE{t) 
for Ui{t, t'), the time evolution operator for Hi{t), has been used by Lambropolous 
p9[] . Certainly for high intensities E{t) this is very problematic, since one requires 
several terms in the expansion to achieve a reasonable result. 

A more promising approach has been advocated by Perelomov, Popov and 
Terentev g, who took H{t) = Ho + V+z-E{t) = Hi{t) and K{t) = Ho+z-E{t) = 
Ho^i(t). Then H(t) — K{t) = V, thus leading to a power series expansion in 
V. Whenever \E\ > \V\ this seems to be very suggestive approximation and is 
based on the fact that the time evolution operator Uo^i{t,t') for HQi(t) is the 
Gordon- Volkov solution, which is known exactly (see equation after (2.29)). By 
the discussion of section 2, in the Kramers- Henneberger gauge this corresponds 
to a power series expansion in suitable translates of V. 

A combination of these two methods has been proposed earlier in a seminal 
paper on the subject by Keldysh ||^, who took the series for Ui(t, t') with z ■ E{t) 
as a perturbation, but in the second iteration step inserted the time evolution 
operator [/o,i instead of the time evolution operator for Hq + V . In fact, it was 



demonstrated by Davidovich et al. ||3^ that to first order the Keldysh approxi- 



mation and the one of Perelomov et al. precisely coincide. When carrying out the 
same steps in the velocity gauge, i.e. for H{t) = H2{t), one obtains the so-called 
Faisal- Reiss approximation |0. 

We want to point out that all such series expansions are somewhat problematic 
since a proper convergence of the series has not yet been established (the only 
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known case is the Born series in scattering theory at high energies, see e.g. 
and the references given there. For the one dimensional situation convergence 
may be shown for integrable potentials |3^), nor is it straightforward to give 



precise quantum mechanical estimates of the first terms. Most statements seem 
to be based on crude semi-classical estimates [0 or in the belief that features 
which have been observed for relatively simple one dimensional models, which 
are anyway put in question [|ll|], carry over in general |jT3|, |10 . 



For realistic pulses for instance with smooth adiabatic turn on and off, the 
first terms will only give reasonable results when the full power is reached but 
will be poor, if not completely invalid near the turn on and turn off point. 

We now give a new rigorous upper bound on the ionization probability. The 
proof of this bound bypasses the problem of summing the whole perturbation se- 
ries by staying strictly with the Du Hamel formula. In what follows the potential 
V will be supposed to satisfy the conditions given in |^ which are tailored to 
ensure the existence of the time evolution operators for Hi{t){i = 1, 2, 3) given in 
(2.19). In particular such potentials V are Kato small (see e.g. 0), i-e. there 
are a < 1 and b < oo such that for all ip in the domain V{Hq) of Hq = —A/2 

||\/^|| <a||-AV^||+6||V^||. (3.1) 

Also the domain 'D{H) of H coincides with T>{Hq). We note that the potentials of 
atoms or molecules arising from Coulomb pair potentials belong to this wide class. 
Also all potentials (except the (5-potential) like smoothed or shielded Coulomb 
potentials used in numerical computations in this context are Kato small. Also 
Hamiltonians with Kato bounded potentials (3.1) are bounded below and if the 
(pair) potentials decay suitably at infinity, then there are no positive eigenvalues. 
Indeed, it has been shown ^ for a large class of potentials including atomic 
and molecular ones that the eigenvalues are contained in [inf a{H),0] where 
cr(if ) denotes the spectrum of H. With these specifications on V in mind we are 
now in the position to state, prove and comment on the first main result of this 
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section. 



Upper bound 1: Let Tp he a normalized hound state of H = Hq + V with 
energy E < 0. Then for any pulse E{t) with l/26(r)^ < ~E, the ionization 
prohahility satisfies the upper hound 

\b{r)\ 



lW-^< / \\iVix-cit)e,)-Vix))mt+\cir)\\\p,iP\\ + 



h{r) 



\Pz 



(3.2) 



Note that the condition on E{t) just says that the classical energy transfer of 
the pulse is less than the classical ionization energy. 



The proof is based on a combination of arguments used in p5| and and 
goes as follows. Since 

exp —itHip = exp —itEip , (3.3) 
by using the Kramers-Henneberger transformation we have 

/(^)^ = ||(l-P)t/i(r,0)^|| = 11(1 -P)T(r)f/3(r, 0)^11 
= 11(1 - P) exp-ib{T)z ■ expic{T)p^ ■ Us^t, 0)ip\\ 
< 11(1 — P) exp ~ih{T)z ■ exp ic{T)pz{U3{T, 0) — exp —iTH)tp\\ 

+ 11(1 — P) exp —ih{T)z ■ expic{T)pz ■ ip\\. (3.4) 

We start with an estimate of the first term on the r.h.s of (3.4). Obviously it is 
bounded by 

||(t/3(r,0)-exp-M)z^||. (3.5) 
We now invoke Du Hamel's formula to rewrite (3.5) as 

r U'i{T,t)[V{x-c{t)e,)-V{x)]exp~i{T-t)H-i)dt . (3.6) 

JO 

Now we use the unitarity of f/3(r, t) (besides the fact that we never iterate 
Du Hamel's formula, this is the crucial step in avoiding perturbation theory) and 
(3.3) to estimate (3.6) by 



\\{y{x-c{t)e,)-V{x 
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dt 



(3.7) 



which is the first term on the r.h.s. of (3.2). To estimate the second term in 
( |3.4|) , we use the triangle inequality to obtain 



— P) exp —ib{T)z ■ expic{T)pz ■ ip\\ < \\{exp{ic{T)pz) — 1 ■ ip\\ (3.8) 



valid for any selfadjoint operator to estimate the first term on the r.h.s. of (3.8) 



which is the second term on the r.h.s. of (3.2). It remains to estimate the second 
term in (3.8). By assumption {1—P)H > 0. Hence for any 5 > 0,{l—P){H+6)^^ 
exists and is norm bounded by 1/6. Therefore 

P)exp-ib{T)z-ip\\ 
= - P){H + 6)-\H + 6) exp -tb{T)z ■ -^W 

= P){H + 5y^exp-ib{T)z-exp ib{T)z ■ (H + 5) ■ exp -ib{T)z ■ 

< hexp ib{T)z ■ {H + 6) ■exp-ib{r)z -^l (3.10) 



We now use the fact that 



+ P)exp-ib{T)z -tpW 



Now we use the estimate 



||(expzA-l)^||<||A^|| 



by 



(3.9) 



exp ib{T)z ■ H ■ exp —ib{T)z 



\{-^V-b{r)e-,f + V 
H - b{r)p, + h{rf . 



(3.11) 



Inserting this into (3.10) gives 



II (1 - P) exp -tb{T)z . < ^-\\(E - b{T)p^ + hirf + 6)^\\ 



(3.12) 



Making the choice 



S = -E-- 




(3.13) 
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which by assumption on 6(r) is > and inserting into (3.12) gives the third term 
in (3.2) concluding the proof of the upper bound. 
We now comment on this result. 

Inspection of the proof of the main result in shows that one has the 
alternative 
Upper bound 2: 

liij)-^ < r\\(V{x-c{t)e,)-V{xM\dt+\c{T)\ ||p,^|| + |6(r)| \\z^\\ (3.14) 
Jo 

which differs from (3.2) only in the last term. Typically near threshold, i.e. for 
small \E\, both \\zip\\ and l/{—E — ^^(t)^) become large, whereas 

WpM = < (2(z^,i/o^))^ 

stays finite. Thus for V being the Coulomb potential 

{i:,Hoi:) = -E (3.15) 

by the virial theorem (see e.g. |^2[). For s-states ipnoo one can improve slightly 
since ||p^V^nooP = |(V'noo, ^^oV^noo) = Similarly one has Wzi^W^ < {ip,r^^){= 
^[Sn^ + 1 - 3i{i + 1)] if tfj = ipnEm, see e.g. @). Again for s-states H^^P = 
r^?/)), which is a slight improvement. Therefore (3.2) and (3.14) are essen- 
tially equivalent. 

We now discuss the first two terms in (3.2) and (3.14). In general, for Kato 
bounded potentials V{x){—A + is a bounded operator. Also since —A is 
translation invariant, we have 

\\V{x-y){~A + l)-i = \\V{x){-A + (3.16) 

In particular for the choice of the Coulomb potential, we prove in Appendix A 
that 

^(-A + 1)"^ <6.35. (3.17) 
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Thus for the general potentials V considered, the first term in (3.2) and (3.14) is 
bounded by 

r||\/(^)(-A + l)-il| ||(-A + 1)^11 + r\\Vn = r\\V^\\ (3.18) 

+ r||V(f)(-A + l)-i ||(2i7o + l) 



which involves E{t) only through its duration but not its strength. For the 
eigenfunctions ipnim of the hydrogen atom, one has (see e.g. 



\\{2Ho + l)t/Jnemf = ||((2/7+l)-2\/)^„,^f 

= {2En + 1)^ + 2(2E„ + l){tpnem, -i'nim) + Mi'nem, Xi'nem) 

= 1 7 + — TT- 3.19 

This quantity is < 8 and behaves like l + 0{^) for n large uniformly in < £ < 
n — 1. Hence the right hand side of (3.18) for bound states of the hydrogen atom 
is bounded by 19. 4r uniformly in n and by 6.35r for all large n. 

For the ground state ^/'loo of the hydrogen atom the first term in (3.2) and 
(3.14) has a much better estimate. As shown in Appendix B 

||(y(f-y)- V(f))^ioo|| <2 (3.20) 

holds for all y G IB? such that the first term in (3.2) and (3.14) is now bounded 
by 2r. 



By a theorem of Pfeifer [38] for the survival probability \{'(p,Tpr)\ of a state 
with ipr = U{t, 0)ip for any time dependent Hamiltonian H{t) one has for all 
small r (see for precise conditions) 



V'r)| > cos( / A(t)dt) (3.21) 

JO 

where in the present context with H{t) = Hi(t) (see (2.19)) 



A{t) = \E{t)\-a^ 

= i\\zi;f-{^,z^)'y^<\\zi;\\. (3.22) 

19 



This gives for the ionization probabihty 



m < l-\{^,^,)\'<ali^£\E{t)\dt^ 



2 



2 



< ^J^ \Eit)\dt) Wz^pr- (3.23) 

This may be compared with the discussion above. (3.23) is weaker than (3.2) and 
(3.14) in the sense that it does not show independence of the field strength when 
6(r) and c(r) are small or even zero. Otherwise it is basically equivalent to (3.2) 
and (3.14) or even stronger whenever the last two terms there dominate. We note 
that the rigorous bound ( p.23| ) may be compared with first order perturbation 
theory in E{t) which gives 



41 = ^1^ Eit)dty\\pz^pr< 



E{t)dt 



2 



■ (3.24) 



We now turn to a comparison with other approximation methods based on 
perturbative expansions used in this context. Since all these approximations 
resolve around the same principle, i.e. an expansion involving the Gordon- Volkov 
time evolution operator f/o^i(t,t'), we will mainly concentrate on a recent work 



by Geltman [13|, who presented an explicit, partly analytical, partly numerical 
analysis of the full three-dimensional hydrogen atom. Also a nice discussion of 
other works may be found there. Geltman employs the approximation method 
of Perelomov et al. in order to compute the excitation and ionization rates for 
the ls/2s/3s/2p/3p/3(i states of the hydrogen atom struck by a linearly polarized 
monochromatic laser pulse of the form E{t) = Eq cos cut. The value for the electric 
field strength in atomic units is chosen to be Eq = 5, 10, 20 and the frequency 
u = 1.5 . Geltman obtained the following general features 

a) At integer cycles, that is for r = 27m /u the rate of ionization becomes 
independent of the electric field strength. In particular for the non-s-states 
it goes to zero. 
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This is reflected qualitatively in our results in the following way. At integer cycles 
6(r) = c(r) = such that the last two terms in (3.2) and (3.14) vanish. Also 
the first term is independent of the field strength, however, not zero and by the 
above discussion too large for the above choices of Eq and u. 

b) The maxima of 6(r) and c(r) are located at half-integer cycles, i.e. r = 
27r(n + i)/cj. 

For the applied pulse the bounds (3.2) and (3.14) also reproduce this feature 
qualitatively but again we emphasize that these bounds hold in more generality 
for all Kato potentials and all states. 

Stabilization for strong, short electric fields has been a highly controversial 
issue with disagreeing results between numerous authors on one side, for a review 
see ||20|, as well as on the other side, see [|19| for a review on these. 

The following result shows absence of stabilization for sufficiently strong, short 
pulses, namely when 6(r) becomes large and r small the ionization probability is 
close to 1. 

Lower bound: Let ip be a normalized bound state of H = Hq + V with energy 
E < 0. Then for any pulse E{t) with l/26(r)^ > —E the ionization probability 
satisfies a lower bound of the form 

Ii^)>l - \ r \\iVix-cit)e,)-V{x))ilj\\dt (3.25) 







■UVix - c{r)e,) - V{xm + ^ \\pM 



Note that now the condition on E{t) is that the classical energy transfer of the 
pulse is larger than the classical ionization energy. Recall that by our previous 
discussion the norms appearing in the first two terms in the bracket may be 



estimated independently of the field strength, such that ( |3.25|) gives a bound 
which involves E{t) only through r and &(t). 
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We turn to a proof. In order to obtain a lower bound on 
/(V^) = 11(1 - P)f/i(r, 0)^||2 = 1 - \\PU^{t, 
it suffices to obtain an upper bound on ||P?7i(r, 0)'?/'||. First we write 



||Pf/i(r, 0)^11 = \\PT{t)Us{t, 

= \\P exp —ib{T)z ■ expzc(r)p2 ■ f/3(r, O)'0|| 
< ||Pexp — ^^(r)^ ■ expic{T)pz ■ (f/3(r, 0) — exp —iTH)ip\\ 
+ ||Pexp —ih{T)z ■ exp ic{T)pz ■ ^||- (3.26) 

The first term on the r.h.s. is estimated by (3.5) yielding the first term in the 
bracket in ( ^ ). 

The second term in ( |3.26| ) is treated as follows. By assumption PH < 0. Let 
5 > be arbitrary. Then P[H — 5)^^ is a well defined operator with operator 
norm < 1/5. Hence 

||Pexp —ih{T)z ■ expic{T)pz ■ ipW 
= \\P{H - 6)-\H - 6) exp -ib{T)z ■ expic(r)p, ■ ^\\ 
< \\\{H -6)exp-ib{r)z-expic{T)pz-'4j\\. (3.27) 



In analogy to (3.11) we now use the relation 

exp —ic{T)pz ■ exp ih{T)z ■ H ■ exp —ih{T)z ■ exp ic{T)pz 
= ^(-^V - h{T)e,f + V{x - c(r)e,) 

= H- h{r)pz + \h{Tf + V{x - c(r)e,) - V{x). (3.28) 



Inserting this into (|3.271) we obtain 



\P exp —ih{T)z ■ expic{T)pzip\\ < -||(V^(x — c(r)e2) — V"(a^) 





+ \\\{E-h{r)pz + \h{rf-6)n (3.29) 
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We now make the choice 



5 = E + h{Tf (3.30) 



which by assumption on 6(r) is > and when inserted into ( |3.29D immediately 



yields the remaining two terms in the bracket of ( p.25 ), thus concluding the proof 
of the lower bound. 



on 



We note that these two theorems are compatible with the result in |25 
the Stark kick, i.e. E{t) = Fo5(t). There it was shown that for fixed ionization 
probability of any bound state ipnim of the hydrogen atom Fq scales like 1/n, as 
predicted by Rheinhold et al. M 



We now return to a comparison between our results and those obtained by 
employing approximation methods based on perturbative expansions and we will 
include the lower bound in the discussion. We stress once more the point that the 
lower bound definitely excludes the possibility of stabilization of the bound states 
for increasing electric field strength, when the applied pulse is short in duration, 
since (|37^), despite the dependence on c{t), may be estimated by a constant, say 
C, independent of the electric field. Hence 

lim Utfj) > . (3.31) 

This shows clearly that the electric field has no stabilizing effect and we are 



therefore in disagreement with [12, 15, 14 1. The lower bound also reproduces 



the result obtained through an expansion around the Gordon- Volkov solution, 
namely for monochromatic linearly polarized laser pulses at integer cycles, i.e. 
6(r) = c(r) = 0, the ionization probability becomes independent of the electric 
field strength. We may qualitatively relate the term proportional to to a term 
also observed in perturbative expansion methods and which is interpreted as the 
spreading of the wave. 

To illustrate our results further, we consider now the concrete example of the 
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hydrogen atom. For the ■0ioo state we obtain as our best estimate 

/Wioo) < (27 + \b(T)\ + ^\c{r)\] (3.32) 

/(« > l-(2r + ^ + -|_i^)\ (3.33) 

Taking the pulse to be of the form Ez{t) — Eq cos cot for < f < r and zero 
otherwise, we have 

\b{r)\^^\smuT\ and \c{t)\ = sin^ (^^^ (3.34) 



■ i ■ ■ ■ I . 




time 

Figure 1: Upper and lower bounds for the ionization probability for the ■j/'ioo-state 
of the hydrogen atom in the first cycle of an applied field Ez — Eocos{1.5t) for 
Eq = 5, 10, 20, when neglecting the "spreading of the wave". 



Figure 1 shows a plot of one cycle, when neglecting the term which is inde- 
pendent of the electric field strength, i.e. the term 2r, in the upper (3.32) and 
lower (3.33) bound. At the far ends one observes the curves for the upper bounds, 

24 



which approach the vertical for increasing Eq. That is starting from the outside 
and going inwards the sohd hne = Eq = 20, dotted hne = Eq = 10 and the 
next sohd hne = Eq = 5. Next we have the lower bounds for Eq = 20 = dotted 
line, Eq = 10 = solid line and Eq = 5 = dotted line. We have used the same 
values as in ||13| and compare with figure 3 therein. Figure 1 clearly reproduces 
the features of Geltman's results, indicating the dips at the half cycles and pro- 
ducing an increasing ionization probability for increasing field strength. We do, 
however, not observe any crossing for different field intensities. When including 
the 2r-term, this pattern will be moved above the trivial bound 1. This may be 
avoided when achieving a better estimate for the factor in front of r, for instance 
when integrating explicitly ( |3.7| ) for a given pulse Figure 2 shows the up- 
per bound for four cycles and reproduces the well known oscillatory behaviour 
superimposed by a spreading of the wave-packet of the Gordon- Volkov solution, 
the so-called over-the-barrier ionization. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

time 



Figure 2: Upper bound for the ionization probability for the z/'ioo-state of the 
hydrogen atom in the first four cycles of an applied field Ez = 10cos(50t). 
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4 Conclusions 



In conclusion we can say that, according to our arguments atoms do not be- 
come resistant to ionization when exposed to short ultra-intense laser pulses. We 
therefore disagree with the opposite point of view, which is sustained through 
numerous numerical studies partly based on explicit solutions of the Schrodinger 
equation and partly based on perturbative methods. In particular we have com- 
mented above on the problems of the latter methods. It is not the intention of 
this paper to discuss the problems of numerical methods, but we like to remark 
that those studies are in general very complex and subject to many possible errors 
which are difficult to check for second parties. We think that the virtue of our 
arguments is that they are analytic and transparent to the reader. An extension 
of our results to multi-particle systems (thus including atoms and molecules with 
several electrons and not necessarily electrically neutral) may be found in . 
Needless to say, since our results are of a qualitative nature, in the sense that 
they merely provide bounds and that, since there are no explicit solutions for 
the Schrodinger equation available, for precise predictions of ionization rates one 
needs more numerical data. 
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Figueira de Morisson Faria for useful discussions and to the International School 
for Advanced Studies (Trieste, Italy) for its kind hospitality. 



Appendix A 

The aim of this appendix is to prove the following bound for the Coulomb po- 
tential V{x) = -i(r = |f|) on L^{R^,d^x) = L^. 
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Lemma A: The following operator norm bound holds 



r 



< 6.35 



(A.l) 



The proof optimizes well known a priori bounds which we take from 



Let L°° be the space of all Lebesgue measurable functions (f on such that 
\ip{x)\ < M < oo almost everywhere. The smallest such M is denoted by Hv^Hoo- 
Then one has the a priori estimate ([Q p. 56) 



llv^lloo < a{p)\\ - AV9||2 + &(P)II'^I|2 

where || II2 denotes the -norm, i.e. \\(p\\2 = J \(p{xy\d^x . Here 

a{p) = cp-^ 

Hp) 



with 



1 



cp 



rdX 



and p > may be chosen arbitrary. 
Let i? > be arbitrary and write 

1 

r 

with 

e{r < R) 



e{r > R) 



/\^h + (h{p) + \\V^^\ 



m\2- 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



Then one has the a priori bound (see [^] p. 165) 
1 

r 2 

Since the estimates \\f\\2 < || (— A + l)v^||2) || — Ay'lh < || (—A + l)y9||2 are trivially 
valid, this gives 

1 



-A + 1) 



<aip)\\Vl'\\2 + bip) + \\V^ 



R\ 



(A.6) 
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Obviously 



W^W - — 



(A.7) 



Inserting (A. 3), (A. 4) and (A.7) into (A. 6) gives 

|i(-A + < nV2p-'R^ + + i (A.8) 

for all p > 0, i? > 0. The claim now follows by optimizing w.r.t. p and R. A 
short calculation finally gives 



< 11 



_7_ 

TTii 



(A.9) 



which is (Al). 



Appendix B 

In this section we will study the quantity {ip\V{x — y)''\ip) ior k — 1 and 2 where 
V{x) — — e -K^) is the Coulomb potential and where t/j is the normalized 
ground state wave function ipiQQ for the hydrogen atom, which is rotationally 
invariant. Therefore this quantity depends on \y\ only. 

Lemma B: Both —{ip\V{x — y)\ip) and {ip\V{x — yy\ip) are decreasing functions 
of\y\ fori/j^ Vioo- 

Intuitively this result is clear: ip{x) has its maximum at x = and —V{x — 
y) >0 has its singularity at x = ^ so their overlap is maximal when y = 0. 

Before we give a proof, we first establish an important consequence. Indeed 
we claim that for ip = ipioQ 

\\iV{x-y)-V{x)M<2 (B.l) 
for all y e IR^. To see this we write 

\\{V{x-y)-V{xM\^ = {^,V{x-y)'^)-2{^,V{x-y)V{x)^P) + {ij,V{xfij) 

< {iP,V{x-yfiP) + {iP,V{xfiP) (B.2) 
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since V{x — y)V{x) > as an operator. By Lemma B, the right hand side of 
(B.2) takes its maximum at ^ = proving the claim since = 2 (see e.g. 

El)- 

To prove the lemma, it suffices to consider y to be of the form y = ce^ with 

c > 0. Now we use the well known formula 

1 1 °° /r \^ 
r = -E - P^i^os^) (B.3) 

where r> = Max{r,c),r^ = Min{r,c) and where (r = \x\,^,(p) are the polar 
coordinates of x. By the orthogonality relations of the Legendre polynomials and 
since ip is the ground state we obtain 

-{^\V{x-ce,m = (^1^1^) 

= 4 / r^e-^'—dr (B.4) 
Jo r> 

4 fc poo 

= - r'^e~'^''dr + 4 / re~'^'^dr. 

C Jo Jc 

This function is differentiable in c for c > and its derivative is easily seen to 
be < 0, proving the ffist claim. Next we have (again by the orthogonality of the 
Legendre polynomials) 

Now for < X < 1 we have 

oo 1 -I oo 2^+1 1 

y x^^—n = -y —r. = — (ln(l + x)- ln(l - x)). (B.6) 

2£ + l x,t^2£+l 2x^ ^ ^ ^ ^ ' 

Inserting (B6) into (B5) gives 

{i)\V{x-ce,f\i)) = ^^%e-2'^ (^In (^1 + -In (^1-0^ cir 

+ re-^r 1^1^ (l + ^) - In (l - ^) ) dr. (B.7) 
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Now the right hand side is not differentiable in c. To remedy this we regularize 
and consider the quantity (0 < £ < 1) 



0<F(e,.) = ^ = ?/%e- 
c c Jo 



ln(l + ^)-ln(l-^(l-£)) 



dr 



+- / re-''"[ln(l + -) - ln(l - -(1 - e))]dr. (B.8) 

C Jc T T 



Since hm^^o -^(c, £) = {^l)\V{x — cezy\ip) it suffices to show that 



(B.9) 



for all < £ < 1, since then (B7) is also monotonically decreasing in c > 0. Now 

1 



re 



-2r 



+2 / re 



-2r 



r 1 r(l — e) 
1 1 



r 1 + 



+ 



i-^d-e) 



dr 

dr. (B.IO) 



The ffist integral on the r.h.s. of (BIO) is negative. By (B8), (B9) and (BIO) it 
therefore suffices to show that for any r > 

c 1 



r 1 + 



ln(l + -) < 0. 



and 



Now 



' . +lu(l-£(l-.)l<0. 

r 1 — -(1 — s) \ r 



(B.ll) 



(B.12) 



r 1 /"^ 1 

ln(l + i) = - / -,dc' > 

^ rJo 1 + ^ ~ 



1 1 



+ £L - r 1 + ^ 

which is (Bll). (B12) is proved in the same fashion. This concludes the proof of 
Lemma B. 
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